Macroscopic Quantum Self- Trapping in Dynamical Tunnelling 
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It is well-known that increasing the nonlinearity due to repulsive atomic interactions in a double- 
well Bose-Einstein condensate suppresses quantum tunnelling between the two sites. Here we find 
analogous behaviour in the dynamical tunnelling of a Bose-Einstein condensate between period-one 
resonances in a single driven potential well. For small nonlinearities we find unhindered tunnelling 
between the resonances, but with an increasing period as compared to the non-interacting system. 
For nonlinearities above a critical value we generally observe that the tunnelling shuts down. How- 
ever, for certain regimes of modulation parameters we find that dynamical tunnelling re-emerges for 
large enough nonlinearities, an effect not present in spatial double-well tunnelling. We develop a 
two-mode model in good agreement with full numerical simulations over a wide range of parameters, 
which allows the suppression of tunnelling to be attributed to macroscopic quantum self-trapping. 

PACS numbers: 03.75.-b, 03.75.Lm, 05.45. Mt 



The transition from the classical to the quantum world 
is a subject of intense interest. In particular, the topic of 
quantum chaos studies systems which exhibit chaotic dy- 
namics in the classical limit of H —> [H-Q ■ An important 
phenomenon in driven one-dimensional quantum systems 
is dynamical tunnelling, first identified by Heller and 
Davis [1] . This is a classically forbidden process whereby 
particles trapped in a regular region of phase space may 
quantum-mechanically tunnel to another. The behaviour 
of such systems has provided im por tant insights into the 
quantum-classical transition [7-14|: in particular, the 
period of the dynamical tunn elling is strongly affected 

4|. Dynamical tun- 



by a number of subtle effects [11 



of the tunnelling period with U that we find numeri- 
cally. While previous work reported detrimental effects of 
nonlinearities on dynamical tunnelling [iif , a connection 
with MQST was not made. Surprisingly, at higher non- 
linearities with U > U CI it we find some parameter ranges 
where dynamical tunnelling reappears. This effect has no 
analogue in bosonic Josephson junctions, where MQST 



321 and demon- 



nelling has m ostly been studied in the single-particle 



regime 



IH 15 -Hi, and has been demonstrated experi- 
mentally with ultra-cold atoms in modulated optical lat- 
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21| . Recently it has been shown in less classical Hamiltonian 



tice potentials 

[lj) that atomic interactions in trapped Bose-Einstein 
condensates (BECs) can have a detectable effect for ex- 
perimentally realistic parameters. 

Here we investigate the effect of repulsive atomic in- 
teractions on the dynamical tunnelling of a trapped BEC 
through the variation of the nonlinearity, U. A priori, 
the effect of nonlinearity on a given dynamical system 
is not clear. It has been shown to suppress transport 



has been extensively studied 
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strated experimentally. (3fj| |. 

We begin by reviewing the dynamical tunnelling 
of ultra-cold atoms. For classical atoms in a one- 
dimensional (ID) potential to exhibit chaotic dynamics, 
the potential must be both driven and anharmonic. The 
experiments demonstrating dynamical tunnelling used a 
modulated sinusoidal potential provided by an optical 
HED] . Here ad consider the dimension- 
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in the kicked rotor and oscillator [22|-|25j, Landau- Zcncr 
tunnelling in optical lattices [26|, b ut in contrast enhance 
transport in quantum rachets [27] . 

We find that dynamical tunnelling also occurs for the 
interacting system with U > up to a critical interac- 
tion strength U cr i t . Beyond £/ cr it we find that dynamical 
tunnelling mostly ceases. We connect the dynamical tun- 
nelling suppression to the phenomenon of macroscopic 
quantum self-trapping (MQST) using a two-mode model 
based on Floquet tunnelling states. It allows us to pre- 
dict the critical nonlinearity U CI it from knowledge of the 
nonintcracting system, and to understand the increase 



where k is the potential strength, e the amplitude of the 
modulation, and x and p arc position and momentum co- 
ordinates respectively. Potentials V(x, t) as in Eq. ([T]) can 
be realised on atom-chip traps in the radial direction [lij . 
This potential has the conceptual advantage of not being 
periodic in space; however, the physics we describe below 
will be generic for any one-dimensional potential where 
dynamical tunnelling is realised. 

For our quantum treatment of the system we con- 
sider a BEC subjected to this single-particle Hamilto- 
nian. We assume that mean-field theory is valid and the 
BEC is well described by the wave function ip = ip(x) 
that evolves according to the Gross-Pitaevskii equation 
(GPE) [aj: iH e s§t $ =[ H + U W'\ 2 ] V>, with / dx\tp\ 2 =l. 
Herep of Eq. ([T]) becomes p = ~ih c ffJ^ and consequently 
[x,p] = iHeS. U parametrizes the nonlinearity, stemming 
from s-wave interactions, and h c g denotes the effective 




FIG. 1: (color online) (a) Poincare section for k — 2.3, e = 0.3, 
showing islands of regular motion 7+ , 7_ separated by a region 
of chaos, (b) Husimi function Q(x,p) of the even tunnelling 
Floquet state: Q(x,p)[*] = |H*}| 2 /(27rft cff ) for h cS = 0.5, 
where |a) is a coherent state centered on momentum p and 
position x. 



Planck's constant. It arises naturally when rescaling all 
variables in the GPE to be dimensionlcss 
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and indi- 
cates how "quantum" the system is, with S c ff — > being 
the classical limit. 

The classical system ([!} is integrable for e = 0. The 
Kolmogorov-Arnol'd-Moser (KAM) theorem 0] states 
that regular regions of motion persist in phase-space for 
e > 0, but become increasingly destroyed as e is in- 
creased [3] ■ An example is shown in the Poincare section 
of Fig. [T] (a), where co-ordinates of classical motion from 
a large range of initial conditions are plotted stroboscopi- 
cally, i.e. at times t = 2-nn for n G N. A key feature is the 
two large period-one islands of regular motion I± , traced 
by trajectories of atoms moving in phase with the mod- 
ulation of the potential [jof . The KAM theorem forbids 
classical motion connecting these islands. 

Quantum mechanics, however, allows tunnelling to 
occur between these islands. Consider the linear 
Schrodinger equation (U~0) obtained from the quan- 
tized form of the Hamiltonian (JXJ) . To relate the quan- 
tum dynamics of the modulated system to the classi- 
cal phase-space, we use Floquet states denoted \u n ), 
that are invariant up to a phase under time evolution 
through one modulation of period T, and can hence be 
found as eigenvectors of the time evolution operator: 
U(0,T)\u n ) 



exp [-i\ n T/h cff \\u n ) [14 



The operator 

U(0,T) evolves the wave function from time t = to 
t = T. As U is unitary, the quasi- energy X n is real. 
The period-one islands of regular motion occur in the 
Floquet spectrum as a pair of states that are even/odd 
respectively under the transformation p — > —p, and have 
support on both islands, as shown by the phase space 
Husimi function in Fig. [IJb). We will label these linear 
tunnelling states \u c ) (even) and |u Q ) (odd). 

An atomic wavepacket that is initially localized on a 
single period-one island is a superposition of tunnelling 
states: \u±(0)) — |u o (0)) ± i\u o (0)), where \u+) is lo- 
cated on the island with p > 0. Using the time evo- 
lution of Floquet states, we have U(0,nT)\u±(Q)) ~ 
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FIG. 2: (color online) (a-c) Stroboscopically sampled momen- 
tum space density demonstrating the transition from (a) dy- 
namical tunnelling to (c) MQST with k — 2.3, e = 0.3 as for 
Fig.[T]and U = {0.7,2.7,4.0} x 10^ 3 respectively, (d-f) Com- 
parison of GPE Floquet state populations (lines, |d±| 2 ) with 
two- mode model results (crosses, |c±| 2 ). Black-solid: |d+| 2 , 
red-dashed: |d-| 2 , black-dotted: |d+| 2 + |d_| 2 , black crosses 
|c+| 2 , red crosses: |c_| 2 . 
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gives rise to quantum tunnelling. Its experimental sig- 
nature is a classically forbidden periodic reversal of the 
stroboscopically sampled atomic momentum as observed 
in [20|, . The quasi-energy splitting of the odd and 
even tunnelling states determines the linear period of 
dynamical tunnelling: Xii n = 27rft G ff/|A c — A Q |. 

For U > the problem is nonlinear, and we cannot con- 
struct the operator U from the evolution of a set of basis 
states. Instead, we find nonlinear Floquet states 31| that 
arc solutions d>„ of 
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(2) 



periodic in the time dimension: <fi n (x,t) = cj) n (x,t + 2ir) 
and vanishing for x —> ±00. A state 4> n (x, 0) will reform 
after one driving period of evolution with the GPE up 
to a phase —E n T/h e s, analogous to the linear case 34 1. 



We only consider the even (odd) nonlinear Floquet states 
localized on the islands, labelled \(p e ) (\4>o))- 

Using \4> c / ) we simulate dynamical tunnelling with the 
GPE and U > 0. For illustration we choose k = 2.4, 
e = 0.3 and h c ff = 0.5, solving for nonlinear Floquet 
states up to U = 4 X 10~ 3 . We begin simulations in 
the state \^) = |</>+) = \<f) c ) +i\<f) ), and evolve it with 
the GPE for 1100 modulation periods. We sample the 
momentum space wave function of the BEC once every 
driving period T, with the results shown in Fig. [5] 
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Before 100 modulation periods with U = 0.7 x 10 -3 
the momentum distribution is little changed [Fig. [2ja)]. 
However, on longer time-scales we observe a complete 
reversal of the system momentum after about 500 modu- 
lation periods, demonstrating dynamical tunnelling. As 
the nonlinearity is increased to U = 2.7 x 10~ 3 the ef- 
fective tunnelling period Tdt increases [Fig. EJb)]. For 
U > 3 x 10~ 3 complete momentum reversal no longer 
takes place, and the population becomes trapped in 
phase-space [Fig. E{c)]. This phenomenon is analogous 
to the cessation of inter-well tunnelling due to MQST in 
a bosonic Josephson junction pt^ - t^ . 

For the system in Fig. [2] tunnelling remains suppressed 
as U is further increased. However, we have discov- 
ered parameter regimes where trapping occurs for a finite 
range of U, but tunnelling returns for larger U, as shown 
in Fig. [3ja). Dynamical tunnelling can then persist to 
high nonlinearities — an effect that is not seen in the 
bosonic Josephson junction. 

To understand these results, we derive a two-mode 
model based on the nonlinear Floquet states. We assume 
that the time-dependent solution of the GPE can be ap- 
proximated by ip(x,t) = c + (t)4> + (x,t) + C-(t)cj>-(x,t) or 
equivalently ij)(x,t) = c c (t)(f) c {x,t) + c (t)cj) (x,t). Mem- 
bers of both pairs are orthogonal by symmetry Inserting 
the latter expansion into the GPE, using Eq. ([2]), project- 
ing out the equations of motion for c e (t) and c (t), and 
finally changing basis to c_|_ (t) and c_(i), we obtain: 

d 

ih cS —c± = Ec± + AEc T /2 

+ %>{A 00 }|c T | 2 c± -Om{A co }|c±| 2 c T 

+ [U eo - 5He{A co }/2 - U cc /4 - U 00 /4}\c±\ 2 c ± 

+ [Om{4„}/2 - C/ oc /4 + U 00 /4]\c T \ 2 c T 

+ [C/ e e/4 + U 00 /4 - U eo - mt{A co }/2}cl 4 



0.015 



+ [i3m{A eo }/2 + C/ce/4 - EW4]4c* 



(3) 

- E a and 



Here E = (E e + E a )/2, AE 
the coupling coefficients are defined as: 
U J dx\(j>i{x, t)\ 2 \cj>j{x, t)\ 2 with {i,j} E {e, o} and 
A eo (t) = U J dx(j) 2 (x,t)4>* 2 (x,t). These coefficients are 
periodic in time with period T. Since the tunnelling os- 
cillations usually take place on much longer timescales, 

we replace the coefficients by their average over one pe- 

?T - 



riod, e.g. Uij(t) -» U tj = y J Uij(t)dt. 

To test the model, we extract the amplitudes d±(t) = 
J dxcj>±(x, t)ip(x, t) from the simulations of the GPE, and 
compare them with those obtained using Eq. ([3]). In 
Fig. HJd-f) we find excellent agreement between the two 
mode model and the full numerical solution of the GPE. 

To analyse self-trapping, we consider the population 
imbalance z = N+ — N~ and relative phase ip = 6- — 9+ , 
where c± = y/N±e i6± with N±,6± £ R, following 
Note that N+ + N- = 1. We find 3m{A 0O } is zero for 
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FIG. 3: (color online) (a) Dependence of the tunnelling rate 
T -1 on the nonlinearity U for fi e ff = 0.5, K = 1.3, e = 0.2, 
showing intermittent MQST for U G [1.4, 2.2] x 10 -2 . Black- 
solid: GPE solution, crosses: two-mode-model [Eq. ©] with 
nonlinear Floquet (black) and linear Floquet (blue) coeffi- 
cients. Red open circles: T^ 1 = \E C (U) — E (U)\/2irh c g. 
The inset shows the same data for the parameters of Fig. [2] 
(b) Dynamical tunnelling at high nonlinearity U = 2 for the 
parameters of the main panel of (a). Black-solid: |d+| 2 and 
black-dotted: \d+ 1 2 + \d- 1 2 from the GPE solution. Red-solid: 
c+| 2 from the numerical integration of the two- mode model. 



all simulations and henceforth neglect this quantity. The 
equations of motion for z and <p from Eq. ([3]) are then the 
Euler-Lagrange equations of the effective Hamiltonian 



Hcs = : |(l-C)+«C l cos(^)+/3Ccos(2(^), (4) 



where C = 1 - z 2 _ and A = (U cc + U so )/4 + mz{A_ co }/2 - 
U co , a = (U ee -U 00 )/2-AE, f3 - U co /2-{U cc + U OQ )/8 + 
«He{A eo }/4. For a = 1 and /3 = Eq. g]) simplifies to 
the Hamiltonian of Ref. [2!| , which analysed MQST for 
a BEC in a spatial double-well potential. 

Following (29| we can find Hamiltonian parameters for 
which dynamical tunnelling cannot occur. Starting from 
2(0) = 1, energy conservation requires that for z(t) = 
to occur there must exist a solution to 
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A = a cos [ip{t)] + /3 cos [2<p(t)] . 



(5) 



The atoms are self-trapped when this equation cannot 
be fulfilled for any ip(t). If we assume |Ai?| 3> \/3\ and 
A_E| ^> \(U CC — U 00 )/2\, empirically justified in most 
cases, we find that tunnelling is impossible if U > U CI n = 
2|A_E|/|Ao|. Here Ao = A/U is an overlap integral be- 
tween Floquet states that no longer explicitly depends on 
U, but implicitly through the shape of \4>i(x, t)\ 2 . We can 
then estimate the critical nonlinearity for self-trapping 
from the linear Floquet states, as they are generally very 
similar to the nonlinear Floquet states for U < t/ C rit- As 
long as the simple expression C/ crit = 2|Ai?|/|Ao| is justi- 
fied, t/crit decreases rapidly for smaller H e g and larger k 
or e, owing to a decrease in AE 14 1. 

Equation ((SJ) does not always predict self-trapping. For 
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p <C a (usually fulfilled) the self-trapping condition is 
(free + C7oo)/4 - U oc + 3<Ke{A co }/2 



1 < 
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(6) 



Aside from AE in the denominator, all terms in the frac- 
tion on the RHS are proportional to the nonlincarity U. 
For U ^> AE the nonlincarity then cancels out, and the 
condition © depends only on the overlap integrals Uij/U 
and A eo /U. These again are only weakly dependent 
on the nonlincarity U through the shape of \4>i(x, t)\ 2 . 
In particular for parameters where <j) (x,t) and cf> e (x,t) 
have a significant difference in mean interaction energy, 
I fee — fool; we w iU expect to see a reappearance of tun- 
nelling at large U. This occurs for values of K < 2; We 
find for n > 2 that U ee ~ U oa [35| . An example without 
trapping at large U is illustrated in Fig. [3] The reap- 
pearance of tunnelling, a striking difference to the spa- 
tial double- well case, arises because the nonlinearity here 
affects both the self-energy of each tunnelling mode and 
the effective mode coupling. 

The ID nonlinearity U can be related to experimen- 
tal parameters by accounting for the details of the con- 
finement geometry [3]. For k = 2.3, e = 0.3 we find 
fcrit = 0.004 (onset of trapping) for A" = 8 atoms. How- 
ever, tunnelling will occur for U — 2 with k = 1.3, corre- 
sponding to N = 4590 atoms [3(| . These disparate values 
for A" highlight the importance of our results for any ex- 
perimental realisation of dynamical atom-chip tunnelling. 
We note that the large h e s used here would require very 
tight trapping potentials [13, HH . 

For the set of parameters in Fig. [3J the two-mode 
model correctly predicts the tunnelling period at small 
U and the first range of MQST. In fact, for many param- 
eter regimes the simplified two-mode dynamics of Eq. ((3J 
makes accurate predictions when AE, Uij, and A co are 
calculated using the linear tunnelling states, i.e. <p\ — > u\ 
(inset of Fig. [3J. For larger nonlinearities, the model 
correctly predicts tunnelling but with a period that is 
too long for U > 0.05. We find that in this region 
the two-mode approximation begins to break down, with 
\d + \ 2 + |rf_| 2 as low as 0.8 for the full dynamics. 

In summary, we have demonstrated the analogue of 
macroscopic quantum self-trapping in a bosonic Joseph- 
son junction for the dynamical tunnelling of BECs. How- 
ever, we have discovered parameter regimes where MQST 
is lifted for large nonlinearities. We have shown that 
most of these features are reproduced by the dynamics 
of a simple two-mode model. An interesting extension of 
our work would be to consider the quantum many-body 
two-mode model, as considered for double- well tunnelling 
in Refs. 37, 38|, or considering heating effects that can 
result from nonlinearities in the presence of driving [39| . 
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